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In order to classify modified gravity models according to their physical properties, we analyze the 
cosmological linear perturbations for f(R, Q) theories (R being the Ricci scalar and Q, the Gauss- 
Bonnet term) with a minimally coupled perfect fluid. For the scalar type perturbations, we identify 
in general six degrees of freedom. We find that two of these physical modes obey the same dispersion 
relation as the one for a non-relativistic de Broglie wave. This means that spacetime is either highly 
unstable or its fluctuations undergo a scale-dependent super-luminal propagation. Two other modes 
correspond to the degrees of freedom of the perfect fluid, and propagate with the sound speed of 
\^ • ■ such a fluid. The remaining two modes correspond to the entropy and temperature perturbations of 

the perfect fluid, and completely decouple from the other modes for a barotropic equation of state. 
We then provide a concise condition on f(R,G) theories, that both f(R) and R + f(Q) do fulfill, 
to avoid the de Broglie type dispersion relation. For the vector type perturbation, we find that 
the perturbations decay in time. For the tensor type perturbation, the perturbations can be either 
super-luminal or sub-luminal, depending on the model. No-ghost conditions are also obtained for 
CZ2 , each type of perturbation. 
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Many modified gravity theories have been introduced in order to explain the acceleration of the universe PJ. In 
particular, the first working model of inflation — the Starobinsky model [2| — made use of such modifications that 
l/"~) ' should somehow appear in any effective Lagrangian for gravity. These modifications are usually considered to be 
negligible at low energies. Yet, we do not know much about the reference scale for gravity. Indeed, there might be 
^ ■ two fundamental scales defined by the Planck mass and the cosmological constant. At this level of ignorance, people 
have considered models of dark energy that change only the gravity sector. It amounts to substitute general functions 
of the Riemann tensor for the Einstein- Hilbert R term in the action 0-Q ■ 

If we consider such a departure from the single R term to be responsible for dark energy, it should become important 
only at late times. But if its contribution to the action becomes stronger and stronger as time passes by, the full 
Lagrangian has to be considered as the fundamental one with its own particle content. This implies for example that 
f{R) theories indeed possess an extra scalar propagating field Q- 

Along these lines attempts based on the other non-zero Lovelock scalar, the Gauss-Bonnet combination Q, also 
appeared. In the so called Gauss-Bonnet theories, the Lagrangian for gravity reads C = R + f(G) 01 ■ Recently a 
paper [PJ showed the presence of classical matter instabilities in such models, ruling them out as a sensible theory for 
gravity and dark energy. Yet, in addition to f(R) and R + f(G) theories, one still has the freedom of introducing 
the most general Lovelock modification of gravity, that is C — f(R,Q). Naively, one may expect that the stable 
perturbation behavior of f{R) can overcome the unstable behavior of f(G)- However, the story is not so simple. For 
the vacuum case (i.e., absence of matter), it was shown in [10j that if the theory gives a non- vanishing JrhIgg ~ /L, 



O 



* 



I. INTRODUCTION 



'Electronic address: dcfclicc@rs.kagu.tus.ac.jp 
^Electronic address: jcan-niarc.gcrard@uclouvain.be 
•f Electronic address: teruaki.suyama@uclouvain.be 



then the nature of scalar perturbations qualitatively differs from the one for a class of theories to which f(R) and 
R + f(G) belong. The study was extended to a system where a single scalar field is (not necessarily minimally) coupled 
to gravity and similar results were obtained [111 ]. 

It is necessary to develop theoretical and experimental tools to disentangle all these alternative theories of gravity. 
We believe they should be classified according to some physical properties. In this paper we argue that the high-fc 
modes behavior, i.e. their scale dependent speed of propagation, in the presence of a perfect fluid (not only in vacuum 
as in [10() can provide such a key physical feature according to which the various theories of gravity can be classified. 

The scope of this paper is to develop the cosmological linear perturbation theory [12[ for the general f(R, Q) theories 
in the presence of a perfect fluid. In particular we will find conditions in order to remove ghosts [13| — by imposing 
the kinetic operator for the fields to be positive definite — and Laplacian instabilities sourced by a negative squared 
speed for the propagating fields. Clearly, this study is relevant since our universe contains matter components which 
are well described by a perfect fluid. Therefore, the results obtained in this paper not only reveal the theoretical 
structure of the model but are also useful for the purpose of constraining modified gravity theories from observations 
of the cosmic structures. 

In section |n] we introduce the action and write down the background equations of motion. Linear perturbation 
theory for the scalars are discussed in Section Mil Vector and tensor modes are discussed in Sections IIVI and [V] 
respectively. Our conclusions are given in Section [VII 



II. THE ACTION 



The action we study is given by 

S = 



16ttG ^' ' 



(1) 



where G is the Newton constant and p is the pressure of a perfect fluid characterized by a chemical potential (/i) and 
an entropy per particle (s). The four- velocity of the perfect fluid is given by potentials, 

u v = -(d v t + 6d v s + Ad v B), (2) 

M 

where £,6, A and B are all scalar quantities. The action (jTJ) for a perfect fluid has been introduced in [1J], and 
its application to the cosm olog ical perturbations for a single perfect fluid within General Relativity(GR) has been 



developed in a recent paper 15| . The fundamental fields that should be variated to derive the equations of motion are 
g^w,£, 0, A, B and s (p is written in terms of these fundamental fields by using the normalization condition u^u^ = — 1). 
The variation with respect to g^ yields the gravitational field equations: 

R^v ~ %9p u ~ M" = O 71 "^ 2/iI/j (3) 

where E^ is the effective energy momentum tensor defined by 

£„„ = V M V„F - g^UF + 2i?V M V„£ - 2g^RU^ - 4R M A V A V„£ - 4i?„ A V A V^ + 4R M1/ D£ 

+ Ag^R afi V a V^ + m mPu V a V^ - \ g^V + (1 - F) (R^ - \g^R) . (4) 

InEq. Q, V = FR + £g-f(R,g), and 

F =dR> t= W (5) 

Meanwhile, the variations with respect to £, 9, s, A and B, with a help of the first law of thermodynamics dp = 
nd/U — nTds (n being the number density and T, the temperature), yield the following equations of motion: 

d a (ny/=g~u a ) = 0, u a d a s = 0, u a d a 9 - T = 0, u a d a A = 0, u a d a B = 0. (6) 

The first equation represents the conservation of the particle number and the second one, the conservation of the 
entropy. From Eqs. © and the relation p = n/J, —p, we can derive the standard energy- momentum conservation [141 ] : 

V^Tt = V p [(p + p)u»u v + p 8^] = 0. (7) 



Before going to the linear perturbation theory, let us first derive the evolution equations for the Fricdmann-Lcmaitrc- 
Robertson-Walker(FLRW) flat universe whose metric is given by 

ds 2 = -dt 2 + a 2 (t)dx l dx l , (8) 

with u M = (1,0,0,0), the four-velocity and H = a/a, the Hubble parameter. In a FLRW universe, all the physical 
variables for a perfect fluid depend only on the cosmic time t: 

P = p(t), p = p{t), n = n(t), T = T(t), /« = /*(*), » = *(*)■ (9) 

On the other hand, the velocity potentials can depend on the spatial coordinates: 

A = A(x), B = B{x), 0= I dt' T(t') + 9{x), (10) 

where A, B and 9 are arbitrary functions leading to equivalent physical backgrounds consistent with homogeneity and 
isotropy. We will take advantage of this freedom to simplify our perturbation studies. 
In a FLRW flat universe, the gravitational field equations ([3]) reduce to 



3^=i 



±V-3HF-12H 3 £ + 8ivGp , (11) 



F = A£H A -4£ET + FH -8£HH -8irG{p + p) -2FH (12) 

and allow us to remove respectively V and F in all the remaining calculations. The continuity equation p+3H (p+p) = 
is a direct consequence of the previous two equations. Indeed, 

R = 6(2H 2 + H), (13) 

G = 24,H 2 (H 2 + H) = 2AH 2 -, (14) 

a 

with y^gG = 8d(d 3 )/di a total derivative with respect to time, as it should be. 



III. SCALAR MODES 

For the scalar gravity perturbations, we expand the metric as 

ds 2 = -(1 + 2a) dt 2 - 2a{t)d l fidtdx l + a 2 (t) (5 tj + 20% + 2d l d J -f) dx l dx j . (15) 

For the scalar matter perturbations, we introduce S£, 69, 5s, 8 A and SB and consider the simplest solution for A, B 
and 9, namely 

A = B = 9 = 0. (16) 

With this particular choice, the scalar velocity perturbation defined by Suj = p~ 1 diV is related to the potentials 
introduced in Eq. ([2]) as v = St + 9(t) 8s. Therefore we also have that 8T°i = ndiv. 

To study the perturbation dynamics, we have to expand the action at second order in the fields. As it is well 
known, the perturbation variables are not necessarily dynamical fields and we can reduce the second order action by 
eliminating the auxiliary ones. After integrations by parts, one is left with a Lagrangian that contains only gauge 
invariant fields and the action for the perturbations can be written as follows 

S {2) = \ fdtd 3 xa 3 [A ab V a V b - a- 2 Be ab VV a ■ W 6 - a - 4 D 11 (V 2 V 1 f - a - 2 E ab VV a ■ W 6 

- Ce ab V a V b - M 2 b V a V b - a- 2 E 13 V8s ■ W x - C^Ss - C 23 V 2 Ss 

- 2 TO 2 3 V 2 Ss - K8s 2 + n {59 M 5s - 5s59 M + SA8B - 5B5A)] , (17) 
with e ab = —e ba and e\2 = 1. In Eq. (Tl7|) , we named V\ = 4>m, V% = %, and 89m with 

H(8F^im (18) 

F + AH^ 

p + P SF + 4H 2 St 
vm = v H -r ——. , (19) 



n F + 4ff 2 £ 
hi = 5e _T(SF + 4H 2 SQ 
F + 4ff 2 £ 



(20) 



The reduced action (fT7|) whose coefficients are given in appendix [Al contains six fields : V±, V 2 , Ss, 59m, 5 A and 
SB. Among them, 5A and SB completely decouple from the other fields. We can also see that SA and SB do not 
contribute to the four- velocity (see Eq. ([2])) since A = B = at the background level. Therefore, they never affect 
physical quantities and can be neglected in the following. 

While V\ and V 2 have terms that are quadratic in their time derivatives, Ss and 56m have only terms that are 
linear in their time derivatives. By taking a variation with respect to 59 m, we find that Ss is a constant of motion, 
Ss = 5s(x). This fact implies that in the equations of motion for the remaining fields, it should be considered as a 
source. The variation for Ss gives the dynamics of 56m as follows 



2nS9 M - 2K5s - 2m 23 v M + ^ V 2 M - C 13 <^m - C 23 v M = . 

a 1 



(21) 



The equation of motion for V a becomes 



_L__d 



a 3 A ab V b 



Be n 



-V 2 H - Ce ab V b 



D 1 



-<5alV 4 ^ 



1 

772 



E ab + \{B + HB)e, 



Ml b -UC + ZHC)e at 



Vh 



2a- 



■5 al V 2 5s 



26 1 2 
Ss d 
2^ di 



(a 3 C a3 ) 



V 2 H 
m 23 5 a2 5s = 0. 



(22) 



Since there is no coupling between 59m and 4>m or vm , we can first solve the equations of motion for the two latter 
fields. Substituting the obtained solution into Eq. (|21[) . we can eventually determine the time evolution of 58m- For 
a perfect fluid of barotropic equation of state p = p(p), we find that the entropy perturbation Ss decouples from the 
other modes in the action (fTT| (see appendix |Bl case I) such that V\ and V 2 form a closed set of evolution equations. 
For p = p(T) and /rr/qq — j\g = 0, then we have instead K = 0, and there is no quadratic term for the entropy 
perturbation in the action (fTTf (see appendix |Bl case II) . 



A. Dispersion relation 



Dispersion relations are necessary in order to understand the Laplacian instabilities for the perturbations. So let us 
derive them for the short wavelength modes for which the time evolution of the background universe can be neglected. 
We find that both B and Dn factors are proportional to the combination (see appendix |A"| : 



A = fRRfgg - f R g — F R ^g - Fg^R. 



(23) 



The qualitative behavior of the perturbations crucially depends on this quantity. Indeed, let us first assume that 
the modified gravity theory satisfies A^O. In this case, after performing a Fourier transformation, Eq. (|22[) can be 
approximated, for large fc's, as 



A ab V b + ^I^V b 



D n k 4 



SalVl + 



UB + HB) e al 



;V b = 0. 



(24) 



From this equation, we can derive the dispersion relations for four modes. The dispersion relation for the first two 
modes is given by 



B 2 + A 22 D n k 4 
A X1 A 22 -A\ 



(25) 



12 



This is exactly the same relation as the one for the non-relativistic de Broglie wave. For the other two modes, the 
dispersion relation is instead given by 



F 22 F>ix 



(26) 



2 ~ B 2 + A 22 D n a 2 ' 
So, we have two propagation speeds, defined as the group velocity aduj/dk, whose squared expressions are given by 



du>i 
l ~dk 

doj 2 
l ~dk 



256 



H 2 A 



k 2 



i F g H 2 + 16 H^g + F R ) {F + 4HZ) 



P 
P 



dp 
dp 



(27) 
(28) 



The Eq. ((27)) is obviously model dependent and can be used to constrain gravity models. Here we want to point out 
that the first two modes have a fc-dependent c 2 , a feature that does not exist in GR and in modified gravity theories 
such as f(R) and R + f(Q). This feature only appears in the general modified gravity theories f(R 1 Q) with A^O, 
and affects only the scalar perturbations. The dispersion relation for those modes has been derived for the vacuum 
case in [HJ, and its expression exactly coincides with Eq. (f27|). The only difference is the change in the background 
expansion of the universe due to the existence of a perfect fluid. This result was expected because a minimal coupling 
of matter to gravity should not affect the high-fc regime. As discussed in detail in [T(J], Eq. (|2T|) shows that, for 
short wavelength modes, either the propagations become super-luminal if c 2 is positive or the spacctimc becomes very 
unstable if c 2 is negative. Here we have explicitly confirmed that the inclusion of a perfect fluid does not alter the 
nature of this physical property at all, i.e. we still have either a super-luminal propagation or a spacetime instability 
if A 7^ 0. Among these two characteristics, the second property might be inconsistent with a viable cosmology and, 
for this reason, could be used to rule out models with A ^0. The first property, the super luminal propagation, is 
not a problem per se as it does not directly violate causality on the cosmological FLRW background. However it 
indicates that the UV completion of the theory may not be Lorentz invariant [lg|. We conclude that f(R, Q) modified 
gravity theories can indeed be classified according to some physical properties: if A 7^ 0, one propagation speed is 
scale-dependent {c\ ex fc 1 ); if A = 0, the propagation speeds are all scale- independent. So direct observations might 
disentangle them. 

The other two modes having a squared speed of propagation c 2 represent the propagation of scalar perturbations 
arising from the perfect fluid. Because of its minimal coupling, the propagation nature of a perfect fluid is not affected 
by the modification of gravity. The remaining two modes only change on cosmological time scale, i.e. w 2 ex k°, as 
shown in the appendix [Cl 

Let us now consider the subset of theories for which the condition A = is fulfilled. It is obvious that both f(R) 
and R + f(Q) theories belong to this class. As shown in [1(|, however, there is an infinite number of other f(R,Q) 
theories that belong to this class. If A = 0, then the B and Dn terms vanish in Eq. (|24[) . Thus we expect that the 
ultra-violet behavior drastically changes, and the de Broglie type dispersion relation should be modified. To see this 
explicitly, let us solve the resulting evolution equations for large /c's: 

A ab V b + E ab ^V b = 0. (29) 

a 

We find that the modes propagate with speeds given now by 

2 116 iH 2 'i - 64 H 3 £ 2 - 64 i 2 HH - 12 £,H 2 F - 16 FH£ - 16 £ff F + 4 £F - 3 FF 



3 16H 3 i 2 +4£,H 2 F + 4:FH£ + FF 



(30) 



c| « S . (31) 

In the case of f(R) theories, £ = and we recover the well-known result c 2 = 1. Theories which are supposed to 
lead to dark energy at late times need to reduce to GR at early times. In this case one typically has fggH 6 <C 1 and 
IrrH 2 <1 0, S|, so that F and £ can be considered as corrections with respect to other GR-like leading terms. In 
this limit, the expression for c 2 becomes 

c? rj 1 + — -. r . 32 

1 3 F + AH 2 £, V 

When this happens, also the background is GR-dominated and 2H ps —3(1 +p/p)H 2 , so that 

c 2 «l- 2 (l + ^)^— . (33) 

If F -C 4iJ 2 £, we are in the Gauss-Bonnet gravity regime, for which evidently the squared speed of propagation 
becomes in general negative. This result, applied to the R + f(Q) theory, confirms the one found in [9|. If F ^> 4i7 2 £, 
then we are in the f(R) regime for which c 2 ~ 1. Therefore, one has to require these theories to be closer to the f(R) 
regime rather than to the Gauss-Bonnet one. 

B. Ghost conditions 

Once we imposed the speed of propagation to be positive, one should also make sure that the modes are not ghost- 
like. This is achieved by demanding the matrix A to be positive definite. For that purpose, one starts with a simple 



field redefinition 

Vi =dWi, 

V 2 =-^CiWi+C 2 W 2 . 
A22 

In this case, the kinetic matrix A for the fields W\$ diagonalizes into 

A = diag(C? Azi dct A, C 2 A 22 ), 

and the additional conditions to impose are 

n 2 p 9216ir 2 G 2 n 2 H i H 2 (p + p) 2 &A 



Ao 



P + P V 
9 



pj 



>0, 



det A = fjHn 2 (F + 4H 2 £) 2 (Fg + 4H 2 &) 2 {F + 4#£) 

where J is defined in appendix [X] There is an extra no-ghost condition 

K > 0, 



A 



(Fg + 4HH g ) 



>o, 



(34) 
(35) 

(36) 

(37) 
(38) 

(39) 



that can be derived once the field Ss is also integrated out (see appendix [C]). It amounts to impose a negative 
coefficient for 5s 2 in the action (1171) . 



IV. VECTOR MODES 



The vector modes treated from an action point of view were discussed for General Relativity in |l5| . Along the 
same lines, we find here the action for the vector modes in f(R, Q) theories. As for the metric, we introduce two 
vector fields such that got = —aGi and gij = a 2 (C'i.j + Cj^). We define then the divergence part of the perturbed 
3- velocity, 5ui of the perfect fluid as u( . In terms of the gauge invariant Vi = —{Gi + aCj), and fixing the gauge 
SB = 0, we find the following action 



S = i / dtdx 3 a 3 



F " l " 4 ^«mMM ,,/ , \ Vn P + P v v , n P + P v,r 
m7rGa 2 \ d i V V ( d i v i) + 2 (P + P) u l C * ^2~ u l u i + 2 — — u i v * 



The variations of this action with respect to the fields Vi and Cj lead respectively to 

- 2 16irG{p + p) v 

V Vi = — a Uj , 

F + 4H£ 

±[(p + p)a 3 uY}=0. 



(40) 

(41) 
(42) 



This result implies that u( is a decaying mode, and so is Vi, in general. 

Using a Fourier decomposition for all the fields, and integrating out both u( and Vi, one obtains the following 
action 



S = 



»/■ 



dtdr a A QdCi, 



where 



Q 



(p + p)k 2 a 2 {F + AH£) 
(F + AHi)k 2 + 16ttG(p + p)c 



Consequently, in order not to have ghosts for any value of fc, one needs to impose the condition 

F + 4H£ > . 



(43) 



(44) 



(45) 



V. TENSOR MODES 

These modes are the easiest ones to study, as no contributions arise from the new scalars degrees of freedom and 
the perfect fluid matter Lagrangian. In this case the action for the iJ^ modes, with H^j = = H l i: can be written 
after decomposing them into the two polarization modes e®,e§. Calling A^ and A§ their corresponding amplitudes, 
one finds 



2 



S=h dtdx 



„3 3 



F + ^Hj . F + 4£_~ 2 

16irG A ® " 16^G^ (VA0J 



(46) 



and a similar action for Aq. Therefore the no-ghost condition coincides with the one given by the vector modes. 
However, these modes do propagate and their speed of propagation reduces to 



2 F + ^ 



4 = — > . (47) 

This is an independent condition, which implies the classical stability of the spin-2 metric perturbations. 

VI. CONCLUSIONS 

We have developed the linear cosmological perturbations for f(R, Q) theories where a perfect fluid is minimally 
coupled to gravity. 

For the scalar type perturbation, we found that there are in general six degrees of freedom. Two modes out of 
six, which do not exist in GR and appear only in modified gravity theories, represent the gravity perturbations. If 
IrrIgg ¥" Irg 1 these modes obey a dispersion relation which is the same as the one for the non-relativistic de Broglie 
wave: to 2 = L(t)k 4 . Correspondingly, the perturbations either propagate with super-luminal speeds if L{t) is positive 
or are highly unstable if L{t) is negative on small scales. The apparent expression of L(t) exactly coincides with the 
vacuum case given in [10|. This result is sensible as a minimally coupled perfect fluid cannot affect the nature of 
the de Broglie type dispersion relation which solely comes from the modification of gravity. The other two modes 
represent density and velocity perturbations of a perfect fluid, which propagate with a sound velocity c 2 s = p/ p, just 
as in GR. The remaining two modes represent the entropy and temperature perturbations of a perfect fluid. These 
modes completely decouple from the other modes for a perfect fluid with a barotropic equation of state. 

For the vector type perturbation, we found that the perturbations decay in time, which is the same as in GR. For 
the tensor type perturbation, the perturbations can be either super-luminal or sub-luminal, depending on the model. 
We found that both perturbations must fulfill the same no-ghost condition. 

The advantage of the action approach presented here is that we can determine the sign of the kinetic terms for the 
perturbation variables. The kinetic terms must be positive to avoid ghost, which put another constraint on f(R,Q) 
theories. Wc found all the no-ghost conditions for the modes. These conditions must be satisfied as well as the 
stability conditions for the perturbation squared speeds in order to construct viable f{R 1 Q) theories. 
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If we define 



Appendix A: Coefficients in the action (|17[l 

T 1 =F + 4JH£ (AI) 

r 2 = F + 12H 2 £ + 2FH (A2) 

r 3 =p(F + 4H 2 £) + 8ttG( P + pf (A3) 

r 4 = F g + 4H 2 & , (A4) 

T 5 = (p+p)f-pT, (A5) 

T 7 = F + AH' 2 i . (A6) 



and 



J = I 6ttG TJ p H A i 2 g + 8nGTl H 2 p F g £ g + ttGT 2 . p ( A + F%) 
+ 768ttGH 3 H 2 [AttG{ P + p) 2 - HpTx] & A, (A7) 

the elements for the kinetic matrix A are 

576H 4 H 2 T 2 r 2 £gA 



A 3Fl 

All -\^GY% 



2 IQ^GHY^p + pY 
1 6 2 



pTiJ 



and 



_ 12nHT 1 

~ — f 2 — 

1 2 



4nG(p + pf 1 23MnGnH*H 2 T 1 T& g A 

- + /il i + 1 6 



V 



pTiJ 



The interaction coefficients between the two propagating fields <f> and v with 8s are the following ones 



The coefficient for 8s is given by 

r 2/-.2rr2 7V2^2T^2i 



K = n T (^\ +n(^] ^^G 2 R 2 H 2 n 2 Y\j g A 



dpj s KdsJ^ P J 

For the Laplacian matrices E and M, it is simple to write down the elements 



(A8) 



&nH(jp + p)T 1 2304^G , n(p + p)g 4 g 2 r 1 r 3 e g A 

Al2 - ^ — + WD (A9) 

A22 = ^P_ + ^ 2 GWH 2 {p + pn G A ^ 
p + p P pJ 



B _ 384irGnH 2 H 2 (p + p)r 6 ij g A 
D 11 = -^%^, (A12) 



(AI3) 



f287rGr 5 nffff 2 r 6 e g A 

^13 = J , ( A14 ) 

4nr!r 5 15367rGnr 5 rir 3 g 3 g 2 g g A 

C-13 = — 1 fr-f; , (A15) 

pi 2 Jp f 2 

2n 2 T 61447r 2 G 2 n 2 (p + p)ff 2 ff 3 r 5 &A 

p pJ 

8TrGn 2 T 5 3072Tr 2 G 2 H 3 H 2 n 2 r 3 r 5 £ g A ., w . 

™ 23 = ^ r-pp ■ (A17) 

pl2 P J 12 



(A18) 



j2 

E 22 = - , Mj' 1 = 0. (A19) 

p + p 

The other elements (i.e., £n, E12, M 2 9 , and M| 2 ), because of their complexity, cannot fit into a standard paper. We 
provide a file to download them 17 1. 



Appendix B: Useful relations for perfect fluids on a FLRW background 

For a perfect fluid we have the first principle of thermodynamics — together with two equations of state n = n(p, s) 
and T = T(p, s) — written as 

dp = ndp — nTds . (Bl) 

Consequently, if p = p(p, s) we have the following useful relations: 



and, by analyticity, 



In a FLRW universe, s = implies that 



By differentiating the identity 



dn\ fdnT\ 



(B3) 



fi = p/n . (B4) 



p = [in - p, (B5) 



and using Eq. (|B4[) together with the conservation of the particle number (n = N/a 3 ), one then recovers the continuity 
equation 

p=-3H(p + p). (B6) 

These relations are quite standard. Let us consider now four limiting cases in f(R.Q) theories. 

Case I: L 5 = 

The combination F5 = (p + p)T — pT appears many times in the coefficients of the action (|I7p . In particular, Ss is 
no longer coupled to V\ when F5 vanishes. So, it is instructive to understand for which perfect fluids it vanishes, i.e., 

(P + P)(J^) p = Tp, (B7) 

since T = T(p, s) with s = on the FLRW background. The Eqs. (|B4|) and (|B5|) imply then that 

"(|).- T - < B8 » 

namely 

T = f(s)p. (B9) 

For a system with two thermodynamical degrees of freedom, this condition is equivalent to having a barotropic 
equation of state. Indeed, if p = p(p, s) and p = p(p, s), we have 

£).--(!).(£),• 

From Eqs. (jBIOj) . (|B5J) . (|B2|) and (|B3|l we respectively infer 



0p\ 



fj. _ v /n 



§f, M ($)„-(&)„ pm..+nT 



i)s 



ft). T -"(«). 



(*). -+"(e).-(«). "(e). "(e 



nT 
. (BII) 



/ S 
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Using now Eq. (|B9[) . we eventually find 



(e)r-nr/ M ; 



=T, (B12) 



*^~ Kfc). 

so that (dp/ds)p = 0, or, equivalcntly, p = p(p). This last statement is a consequence of Eq. (|B1[) , when we change 
the two independent variables from (/i, s) to (/?, s). After doing this, one should consider (i = p{p, s), and the first law 
of thermodynamics implies (dp/ds) p — n[(dfi/ds) p — T\. 

From Eq. (|A16[) . we obtain C23 = — 2n/(s) such that the conservation of the particle number N = a 3 n together 
with the conservation of the entropy per particle s imply dt(a 3 C23) = 0. Consequently, the field Ss decouples from 
both Vi and V2 in the action (fT7j) . Indeed, C13 defined in Eq. (|A15|I and 77123 defined in Eq. (|A17[) also vanish if 

r, = 0. 



Case II: K = 

It is also useful to understand for which equations of state the coefficient K for Ss 2 vanishes, at least in the A = 
consistent case. This happens for 

or, equivalently, 

T—Wf-(%) CBM) 



(f) 



if, again, the system has two degrees of freedom, i.e., if T = T([i, s). The solution 

H = ii(T t s) = Ts + g(T) (B15) 

implies that 

dp = ndfi - nTds = n[s + g'(T)} dT . (B16) 

In order that dp be an exact differential we need then to impose that 

n = n(T,s)= KT ) (B17) 

so that we finally find {dp/ 8s)t = 0, or, equivalently, p = p(T). 

Case III: K = 0, T 5 = 

There is an equation of state for which both K and T^ vanish, even when A ^ 0. This case corresponds to 
/(s) = —r- and g(T) = cT (c is a constant). In this case we find that 

T=-£-, (B18) 

s + c 

and p = p(T) since p = p(p) and p = p(T). A radiation fluid fulfills these properties. 

Case IV: ideal gas with A = 
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It is also interesting to derive K and T5 for a classical monoatomic ideal gas. The thermodynamical quantities for 
this case are given by 

n( M , *) = ^^ (/x - m) 3/2 cxp(| - a) , (B19) 

T(n,s) = l(u-m), (B20) 

pO*,*)=nT, (B21) 

p((M,s) =mn+^nT, (B22) 

where m is the mass of an atom. We have used the Sackur-Tetrode equation for a monoatomic atom. The energy 
density p is calculated from the identity pn = p + p. Assuming a model satisfying A = 0, we find that K and T5 are 
written as 

K=\nT, T 5 = -2mHp. (B23) 

Since T > 0, we have K > 0, so that no ghost appears for an ideal gas, as expected from our real experience. 

Appendix C: Eliminating the auxiliary field 

Let us integrate out Ss from the action. In that case, 89m gets a kinetic term quadratic in its time derivative. It is 
then useful to define a new variable 

V 3 = S9 M -%2-V a , (CI) 

2n 

which can be considered as a redefinition for 89m- Keeping only the relevant higher order terms in k, the action (|17[) 
becomes 

S=\ J didW (A ab V a V b - a- 2 F ab VV a ■ VV b - a - 4 D ab V 2 V a V 2 V b - a - 2 L ab VV a • W 6 ) , (C2) 

where now both a and b run from 1 to 3, and A ab , F ab , D ab and L ab are all 3 x 3 matrices given by 

A n A 12 0\ / B F\ ( Q 

A = | A 12 A 22 , F = \ B , D=\ ) , (C3) 

0A 33 / \ -F / \ 000 



£11 + ^(C 13 + 3HC 13 ) E 12 + ^{C 23 + 3HC 23 - 2m 23 ) j^[E 13 (K + 2HK) - KE 1 



■'• 



L= ( E 12 + ^(C 23 + 3HC 23 -2m 23 ) E 22 | . (C4) 

jfj [E 13 (K + 2HK) - KE 13 ] 



Here 



i33 



n2 f - nEl3 n-n E '* 3 (a\ 



and all the other quantities have already been defined in the appendix A. 

The corresponding dispersion relation is obtained from an eigenvalue equation 

/ k 2 k 4 k 2 \ 

det \-lo 2 A + ilo— F+— D + — L 1=0. (C6) 

\ a 1 a 4 a z J 

This is 6-th order algebraic equation for to. Using the above equations, we confirm that two solutions coincide with 
uii given by Eq. ([25|) and the other two solutions coincide with u> 2 given by Eg. ([26]) . We also check that the remaining 
two solutions are independent of k, which means that these modes evolve only on the cosmic time scale. 

The action (|C2[) leads to the extra no-ghost condition A 33 > 0, which implies K > 0. Although classically the 
equation of motion for V 3 depends on the evolution of V\ and V 2 , its kinetic term should be positive in order to have 
a sensible spectrum for the modes at the quantum level. 
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